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Abstract 



We present the leading order differential and total rates for J/ip 
production at LEP. By leading order we mean all terms of the form 
a s [a s \og(Ml/Ml)] n and < +1 log'(z 2 ) \og m {M 2 Z / , (I + m = 2n - 1), in 
the regions z = 2E^/Mz ~ 0(1) an d z <C 1, respectively. In the inter- 
mediate region we interpolate using the available data. This resummation 
eliminates the 0[a s (M^,) /a s (Mz)] ~ 2 theoretical uncertainties in previous 
calculations. The log(z) resummation results in a suppression of the small 
z region due to coherent gluon emission. Comparing the zeroth moment 
with the LEP data we find the value for the effective octet matrix element 
to be (6^( 3 S"i)) = 0.019 GeV 3 . The theoretical uncertainties are substan- 
tially smaller than those from Tevatron extractions. Using this value of the 
octet matrix element we make a prediction for the first moment of the dif- 
ferential rate and find that the resummed differential decay rate is in much 
better agreement with preliminary data than the color singlet result or the 
unresummed color octet prediction. 



I. INTRODUCTION 



The subject of quarkonium production has gained renewed interest due to the fact that 
it can now be understood from first principles. The production rates are calculated in a 
systematic expansion in a s and v, the relative velocity of the heavy quarks in the rest frame 
of the quarkonium bound state. This is accomplished by working in non-relativistic QCD 
(NRQCD) where the expansion in v is implemented by utilizing the scaling properties of non- 
perturbative matrix elements ]I[ . The use of this effective theory has clarified formal issues 
and allowed for better fits to data. In general, NRQCD leads to larger cross sections than its 
historical predecessor, the color singlet model, largely because NRQCD predicts substantial 
contributions to various cross sections from the color octet channel. This happens even 
though the octet matrix elements are suppressed by powers of v, because the color singlet 
channels can themselves be suppressed either kinematically or by powers of a s . For instance, 
ip' production at the Tevatron can now be well fit if one allows for color octet production 

i- 

The field of quarkonium production itself has matured to the point where we would 
like to go beyond order of magnitude accuracy. Indeed, once the formalism is verified 
quantitatively, it can be used as a tool in other areas in strong interaction physics, such as 
heavy ion collisions and the measurement of the spin dependent gluon distribution functions. 

Presently, the values for the octet matrix element (Og( 3 S\)) and a certain linear combina- 
tion of the (Og ( 3 Pq)), and (Og ( 1 S'o)) matrix elements have been extracted at the Tevatron. 
These matrix elements should be universal in nature, in that we should be able to use them 
as input in other processes to make definite predictions. Unfortunately, hadronic uncertain- 
ties in the Tevatron extraction make this difficult. The extraction depends sensitively on the 
choice of gluon distribution function and factorization scale, as well as how one treats initial 
gluon radiation 0. Extractions of (Og ( 3 5i)) are as disparate as 2.1 x 10~ 3 [||], 2.7 x 1CT 3 
HI, 6.6 x 1(T 3 [§, and 14.0 x 10~ 3 §, in GeV 3 . A process involving smaller theoretical 
uncertainties is needed. Perhaps the cleanest setting to extract a value of the octet matrix 
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element is in prompt J / ip production at LEP, because in lepton initiated processes the the- 
oretical errors are bounded by our computational strength rather than higher twist effects. 
Furthermore, at LEP we can measure the energy distribution as well, thus once we've ex- 
tracted the value of the octet matrix element, we can then make predictions for the moments 
of the rate. This provides a strong test for the color octet mechanism. 

Formally there are two leading order contributions in the a s and v expansion, both 
in the singlet channel, which are of order 0(a 2 s v 3 ). There is a contribution from gluon 
radiation in the singlet channel Z — > ipgg that is suppressed by powers of Ml/ 'El @. There 
is also the color singlet charm quark fragmentation process Z — > ipcc [HU, which has no 
power suppression and thus dominates over non-fragmentation processes, for large E^. Light 
quark octet fragmentation (in which the mother parton does not combine to form part of the 
bound state) is naively of order a 2 s v 7 , down by v 4 ~ 1/10 compared to charm fragmentation. 
However, as it turns out, this channel is enhanced due to the presence of large logs as well 
as a numerical factor of five due to the number of possible quarks that initiate the process. 
Indeed, previous calculations of the J/ip production rate at LEP are dominated by these 
light quark fragmentation contributions [|10],|rTJ. They give cross sections that are of the 
correct order of magnitude when values of the non-perturbative matrix elements are taken 
from the Tevatron fits ||. However, the same logs that enhance the octet channel also put 
the convergence of the perturbative expansion into question. 

The tree-level calculation of the differential cross section in the color octet production 
channel jlOHI^I is enhanced by a large logarithm, ^(Z — ► ip + X) ~ a 2 s \og(M^/M 2 )/z, 
leading to large double logs in the total rate. Since a s log(M|/M|) ~ 1.5, we should 
treat a s \og(M^/Ml) as order one and resum all powers of the large logarithm. With this 
counting, the octet channel is 0(a®v 7 ), on par with the singlet fragmentation contribution. 
More practically, the tree-level calculation has a factor of two uncertainty associated with 
the scale at which a s is evaluated, since a s (M^)/a s (Mz) ~ 2 (this is just a restatement 
that there is a large logarithm). The resummation of the leading logarithms eliminates 
this uncertainty, so the resummation procedure is essential from both a practical and a 



formal standpoint. We therefore calculate the quarkonium differential production rate at 
LEP taking all terms of the form log m (M|/M?) as leading order. We will see that 
this resummation dramatically changes the differential cross section. However, summing 
the above mentioned logs will only yield the correct leading order differential rate if z is 
sufficiently large. When z is parametrically small, terms of the form a s log(z)/z become just 
as important. Furthermore, these logs will also contribute double logs to the total rate given 
that the lower limit on z is IM^jMz- This second type of log, due to soft gluon emission, 
is resummed using a formalism familiar from discussions of jet multiplicities fl3|]l4l| . Thus, 
we split the calculation into two regimes, z ~ 1 and z < 1. We then interpolate between 
these regimes using the data. 



II. FRAGMENTATION FORMALISM 



The tree-level differential rate for color octet J/ if) production is [[□]] 
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where the rescaled if) energy in the Z rest frame z = 2E^/Mz has a physical range of 



2M^/M Z < z < 1 + M2/M|, 



4M 2 /M§) 1/2 , and to the order we work, is 



twice the charm mass. Performing the integration over z leads to the aforementioned double 
logs. In the fragmentation limit, Eq. (|l|) can be simplified to 
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In this limit, the differential rate can be recast as the sum of quark and gluon fragmentation 
processes, 

dT 



dz 



(Z -> iP(z) qq) = 2 C q (^\ z) * D q (fi 2 , z) + C g (jj, 2 , z) * D g ^\ z). 



(3) 



Here, D q ^(z) and D g ^(z) are the light quark (or anti-quark) fragmentation and gluon 
fragmentation functions respectively. The asterisk denotes convolution with the partonic 
production rates, C * D = f z C(y)D(z/y)dy/y . The ji dependence of the fragmentation 
functions is canceled by that of the coefficient functions, C q and C g . All dependence on 
is contained in the fragmentation functions, while all dependence on Mz is contained in the 
coefficient functions. It is this factorized form that will later allow us to resum the large 
logarithms. 

We choose to define the color octet fragmentation functions according to Collins and 
Soper fl5f , 



D g ^\z) = — - *' ° J dx-e- lP+x -^(0\Gt u (0)al(P + ,0)a4P + ,0)G^(0,x-,0)\0), 
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(/i 2 , z) = J dx-e- iP+x ~/ z -Tr color -Tr Dirac (4) 



7 + (0|Q(0) MP + , 0) a^P + , 0) Q(0, x~, 0)|0) 



Here, G^ v is a gluon field strength tensor with color index b and Lorentz indices + and 
u, Q is a quark field, is a creation operator for a ip meson, and d is the number of 
spacetime dimensions. The fragmentation functions are interpreted as the probabilities for 
a parton with momentum k + to decay into a ip with light cone momentum P + = zk + . We 
have chosen to work in the light cone gauge, where eikonal factors usually written to make 
gauge invariance manifest reduce to the identity. An advantage of this definition over the 



alternative |T6| is its consistency with factorization at any subtraction scale fx. 

For the case of quarkonium fragmentation D q ^ and D g _+^ can be calculated in a sys- 
tematic expansion in a s and v by matching onto NRQCD. Any soft divergences which may 
arise due to the semi-inclusive nature of the process cancel in the matching. Since we are 
working to leading order no such divergences occur, and the matching is trivial in the sense 
that there are no corrections to be calculated in the effective theory. The calculation of 
these fragmentation functions are well documented in the literature pl7| , |T8|| , so here we just 
present our results which are in agreement with these previous calculations. 

In the MS scheme we find 
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where 



(6) 



T a is a color generator, while x an d ^ are two component NRQCD spinors. Inserting these 
fragmentation functions into Eq. (Q) and matching onto the QCD calculation, Eq. @, we 



obtain the coefficient functions C q and C g 



C q (fx 2 ,z) = T(Z^qq)S(l-z) 
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Note that there is no physical distinction between what we call gluon fragmentation and 
what we call quark fragmentation, as it is always possible to shift some finite piece from one 
to the other. However, factorization dictates that we match in such a way as to make the 
Wilson coefficients independent of the long distance physics, i.e. M^. 

Once we choose the scale \x to be O(Mz), there are no large logs in the Wilson coef- 
ficients. They have been shuffled into the fragmentation functions and can be resummed 



via the DGLAP |L9| equations. All the leading logs, of order 0[a s log(Mf /M|)], will reside 
in the quark fragmentation function. The contribution from gluon fragmentation is then 
subleading, contributing only at 0(a 2 ) to the rate. We choose to keep this contribution in 
order to reduce the /i dependence of our result, although we do not claim accuracy to the 
level a 2 . There are other terms of this order, arising from both two- loop running and the 
a 2 corrections to the initial gluon fragmentation function, that have not been included. 

In addition to the color octet contribution we include the contribution from the color 
singlet, which is formally 0(a 2 v 3 ), but as mentioned above is numerically smaller than the 
octet contribution. Again using the Collins-Soper definition for the singlet fragmentation 
function we find 
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z) = ||| (Of ( 3 S0> ^^(16 - 32, + 72z 2 - 32, 3 + 5z 4 ) 
which agrees with [§,|l^ . 



III. z ~ 1 RESUMMATION 

The resummation of the log(M|/M|) is accomplished via the usual renormalization group 

analysis of the fragmentation functions. The evolution equations are given by 
dD q (fi 2 ,z) _ a s (/i 2 ) 



{P g ^ qg * D q (fi 2 ) + Pg_+gq * D g (^)} , (9a) 



where the functions P are the standard splitting functions. We solved these equations 
numerically using m c = 1.48 GeV, \i = M z , a s (M z ) = 0.118, and chose a s (M^) to be 
consistent with one- loop running from Mz- 

The contributions to the differential rate from both the evolved octet (solid line) and 
evolved singlet (dotted line) fragmentation functions are displayed in Fig. 1, in units of MeV. 
We used the Tevatron extraction |J to obtain the normalization of the octet fragmentation 
functions. For comparison, the octet contribution without evolution is shown in dashed 
lines, and the singlet contribution without evolution in dot-dashed lines. Resummation 
of the log(M|/M|) terms greatly enhances both the octet and singlet rates at small z. 
However, as mentioned above we should not trust this result in the small z regime. 

At this point we should point out that there are some values of z ~ 1 where the differential 
rate is not trustworthy. When z approaches within v 2 of one, we begin to probe the hadronic 
structure of the quarkonium state, and the expansion in v breaks down . To correctly 

describe quarkonium production at the edge of phase space we must introduce a structure 
function p2]j23| which resums all the large non-perturbative corrections. Indeed, LEP would 
be an ideal place to study these structure functions if there were more data available. In 
the total cross section however, or any sufficiently smeared version of the differential rate, 
such as the first ten or so moments , the expansion in v is well behaved. 
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FIG. 1. Differential rate / or ^ e octet channel with (solid line) and without (dashed line) 
evolution, and for the singlet channel with (dotted line) and without (dot-dashed line) evolution, as 
a function of z = 2E^/Mz- The octet matrix element has been extracted from the Tevatron (see 
text). 



IV. RESUMMATION FOR SMALL z 



The lower limit on z is z m [ n = 2M^/Mz, so that when z 



we need to include 



a resummation of the log(z) terms. Indeed, the generic term in the decay rate has pieces 
of the form a n s +1 \og\z 2 ) log m (M§/M|), (I + m = 2n - 1). We thus treat log(> 2 ) to be 
of the same order as log(M|/M|), and resum all terms of the above form. This problem 
was first encountered in the calculation of jet multiplicities, where it was noticed that the 
splitting functions are highly singular at small z and need to be resummed. The results of 
these calculations led to predictions for the shape of the hadron multiplicties (under some 
assumptions of quark-hadron duality) which fit the data extremely well [p4j. Indeed, the 
results make the striking prediction that there should be a suppression at small z which sets 
in at higher z than would be expected from just phase space suppression. This suppression 
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FIG. 2. Typical ladder diagrams which lead to the leading logs. To include coherence effects, 
angular ordering is imposed such that the angle between two branching partons is smaller than the 
angle of the previous pair. 

is due to angular ordering of soft gluon emission and is a consequence of gluon coherence 



which is naively missed in canonical ladder resummations [25 |. 

Our calculation differs from those previous calculations of hadron multiplicities in that 
we can actually calculate the hadronic fragmentation function in a systematic fashion, in 
terms of some (in our case effectively one) unknown matrix element. Thus, we are concerned 
with the normalization as well as the shape of the differential decay rate. There are several 



different formalisms for handling the coherent gluon emission problem [TJ|,[14],|26j . We choose 



to follow the formalism developed by Mueller p6| , p7| , and we refer the reader to these papers 
for details. 

We begin by noticing that it is the gluon splitting function which is most singular at 
small z, P g _ gg (z) m 2C^/z. Thus, at leading order as defined above, the branching will all 
come from gluon splitting once the initial quarks mix into a gluon. For the moment let us 
consider pure gluon splitting. All the leading logs come from ladder diagrams of the form 
shown in Fig. 2 (without the initial quark box). The coherence issues can be skirted by 
imposing angular ordering on the gluons such that the angle between two branching partons 
is smaller than the angle of the previous pair |28j. This allows us to rewrite the series via 
the integral equation 

zD 9 (t, z) = 8(1 - z) + ^ f 1 ^ f ^z'D 9 (t', z\ (10) 

7T Jz Z' Jt V 

where Ca = 3 and at the end of the calculation t is taken to be t = M 2 /M 2 z 2 in order 
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to enforce the angular ordering. For now, D g (t, z) has been normalized to 1. Iterating this 
equation leads to 



m=l 
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Taking moments of Eq. (|TT|) gives 
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The above sums have a simple closed form. After reinserting the normalization 
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(n - l) 2 + 8a s C A /ir 

7„ is the resummed diagonal anomalous dimension of the fragmentation function. If we ex- 
pand in a s , then it correctly reproduces the most singular pieces of the previously calculated 



two-loop splitting function j29|. 

At small z, the leading contribution to the quark fragmentation function comes from the 
gluon fragmentation function above, convoluted with P q ^ gq (z) « 2Cp/z. Summing up all 
ladder diagrams of the form shown in Fig. 2 leads to 

C F 
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with moments 
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Note that expanding Eq. (|T6| ) to leading order and convoluting with twice the coefficient 
function, Eq. (f7a| ) (for quark and antiquark fragmentation), recovers the tree- level differential 
rate, Eq. (|2]), in the small z limit. 

We can now use the renormalization group to improve Eq. (|TT|). After running we are 
left with, 



D qn {M 2 z /Ml 
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This equation can now be inverted back to z space leading to the result 



D q (M 2 z ,z) 
where 



- (0 * ( 3 Sl)) | a>iMM|)) _^ 



n 



cxp 



— ( log - 

4a [ 



(19) 



2a" (1 ° g -- 6 



2tt 



2tt 



486 \\ C A a s {Mlf \J C A a s (M 2 ^ 
1 1 

1 



2C f 



2C f 



(20) 
(21) 
(22) 

(23) 



b \y na s (M z ) ^ 7ra,(M|) J ' 
and bo is the coefficient of the one-loop beta function. This result was reached in the saddle 
point approximation where 1 — no is small and thus should not be trusted for z values 
larger than z ~ 0.2. Subleading corrections to this result can be systematically included 
by properly adapting the formalism discussed in [p^ , p0|1 , where they were found to be of 
order y/a~ s at the peakQ. We should note that we expect the relative size of the subleading 
corrections to the total rate to be smaller than those for the differential rate, given that in 
the differential rate the subleading terms are down by log(z 2 ), whereas in the total rate they 
are down by log(M|/M 2 ). 



1 We have checked using the results of [27], that the numerical value of the corrections to the 
coefficient function are actually quite a bit smaller than this. 
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V. EXTRACTION OF THE MATRIX ELEMENT 



We now have leading log expressions for the total and differential rates in the small and 
large z regions. Before we interpolate between the two we need to determine their respective 
regions of validity. Let us investigate the size of the contributions we are neglecting in the 
large z region. The generic term in the decay rate has the form log^ 2 ) log" 1 (M|/M 2 ), 
where I + m — 2n — 1. The first term that is not included in the large z resummation is 
suppressed by log(z)/log(M^/M^). Thus a conservative value for the lower value of z is 
~ 0.5. 

The small z resummation is not valid at large z for two reasons. First, we have used the 
saddle point approximation to compute the inverse Mellon transform, which is only valid 
for n — 1 small. Second, we have neglected less singular terms in computing the small z 
fragmentation result. Therefore, we will trust the small z result only up to the peak, where 
n — 1 pa 0, but not for larger z. We believe this to be a conservative upper bound on z for 
this approximation. 

Thus, to obtain the full differential cross section, we interpolate between the small z ^ 0.2 
and large z ^ 0.5 resummations using the data in this region as our guide. We checked that 
varying the interpolation while still staying within the error bars, changes the result by 
at most 25%. Furthermore, we tested the sensitivity of our predictions to increasing and 
decreasing the lower bound on z in the large z regions as well as increasing the upper bound 
on the small z region. These variations do not appreciably change the total rate or the first 
moment. Note that decreasing upper bound on small z regions does drastically change the 
results. However, taking the maximum value of z to be below z pea k is not a reasonable thing 
to do however, given that the peak is the position where the saddle point approximation 
is trustworthy. In addition, the data for hadron multiplicites fits the resummed predictions 
extremely well near the peak J24 . 



In our final result we also included the non-fragmentation corrections given by the dif- 
ference between Eqs. ([!]) and (Q). These are significant at very small z, and contribute to 
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the total rate even when M^/Mz — > 0. 

Since the data includes feed-down from excited charmonium states, the rate should be 
written in terms of the effective matrix elements || 

(^ (n) ( 3 ^i)) = E (^ (m) ( 3 ^i)) BR(Vi(m) - 1>(n)+X), 

m>n 
m>n 



(Of (n) ( 3 ^i)) = E (Of (m) ( 3 ^i)) BR(^(m) - V(n) + X). (24) 



Saturating the excited states with the ■0' and Xcj, and using values for (C$( 3 Si)>, (Of f^)), 
and (Of J ( 3 Si)) from § gives (6^( 3 Si)) = 0.014 ± 0.002 GeV 3 . This number has at least a 
factor of two theoretical uncertainty. Extracting the analogous color singlet matrix elements 
from the ip and 0' electronic widths Q gives (6f ( 3 Si)) = 1.45 ± 0.10 GeV 3 . 

Combining the singlet and octet fragmentation contributions gives the total differential 
rate. Integrating over -0 energies yields a total branching ratio of 



BR(Z - prompt J/0 + X) = ( 1.47 + 0.47 ( °^ Sl)) , ) x 10" 4 , (25) 

1 0.014 GeV 3 1.45 GeV 3 ' 



compared to a branching ratio of (1.93 (dt^S^/O.OU GeV 3 + 0.68 (6f ( 3 S!))/1.45 GeV 3 ) x 
10" 4 from \M- 



The total branching ratio has been measured to be [31,32 



BR(Z -> prompt J/ip + X) = (1.9 ± 0.7 ± 0.5 ± 0.5) x 10~ 4 OPAL 
BR(Z -> prompt J/0 + X) = (3.0 ± 0.8 ± 0.3 ± 0.15) x 10" 4 ALEPH 
BR(Z -> prompt J/0 + X) = (2.7 ± 1.2) x 10~ 4 L3 

BR(Z -> prompt J/ip + X) = (4.4+ 3 ;^) x 10" 4 DELPHI . (26) 

For OPAL and ALEPH, the uncertainties from left to right are statistical, systematic, and 
model-dependent, while for L3 and DELPHI, they are purely statistical. We use the LEP 
average to extract a value for the effective octet matrix element of 

<6?( 3 £i)) = (0.019 ± 0.005 stat ± 0.010 theo ) GeV 3 , (27) 
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FIG. 3. Differential rate ^ as a function of z = 2E^/Mz vs data. The dashed line is the sum 
of the tree-level octet and singlet results and the solid line is the interpolation between the large and 
small z region octet resummation plus the singlet resummation. 

where the first uncertainty is purely statistical, and the second is theoretical. The theoretical 
uncertainty comes from adding in quadrature roughly 30% contributions from perturbative 
corrections suppressed by a s (M^), higher order matrix elements suppressed by v 2 , and 
subleading logs. 

Since the LEP experiments include feed-down from ip' and \Ji we cannot directly com- 
pare our extraction of (Og ( 3 Si)) with those of (0$ ( 3 «Si)) in HdH , but our value of (Of( 3 Si)) 
is comparable with the central value from ||. While our statistical uncertainties are larger, 
our theoretical uncertainties are under good control. Since this cannot be said of Tevatron 
extractions, where theoretical uncertainties dominate, we believe Eq. (p7|) represents the 
most reliable extraction currently available. 

Fig. 3 shows the complete differential rate using our extracted value for (Of( 3 Si)). The 
solid line is the sum of the resummed singlet fragmentation and the interpolation between 
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the large and small z regions for octet fragmentation. The dashed line is the tree-level 



result and the data is from the ALEPH collaboration |32j with efficiency correction from 
f33|j . Given the large errors, it is difficult to make any definite statements, but it seems the 
resummed rate fits the data better at this time. The effect of evolution is to enhance the 
small- z peak, a promising experimental signature of the octet component. A manifestation 
of this signature is a relatively small first moment, which we find to be 

1 m . ! dr ' Z - pI °7 tJ ^ + X ' , * = 0.30, (28) 

T(Z -> prompt J/ip + X) J dz ' K J 

while the tree-level differential rate, Eq. (|l|), gives ~ 0.5. A very rough estimate of this 
quantity obtained from the data |E| suggests a value of 0.26 ± 0.10. This is in sharp 
contrast to the color singlet prediction. The tree-level color singlet decay rate predicts the 
ratio of the first moment over the zeroth moment to be 0.62. Resummation softens the color 
singlet decay rate, but the ratio is still too large, 0.47. The ratio is independent of the color 
singlet matrix element. Therefore, even if the color singlet rate were arbitrarily increased 
so that the color singlet model fit the experimental branching ratio, the ratio of moments 
would not fit the data. A rigorous extraction of the first moment by the experimental groups 
could provide an extremely clean, quantitative test of the NRQCD approach. 

In conclusion, we have resummed large logarithms in the rate for prompt J / ip production 
at LEP to obtain the leading-order prediction. We predict a branching ratio that is slightly 
smaller than the tree-level prediction [|H],|ri|]. Moreover, we have eliminated a factor of 2 
uncertainty in the tree-level result. Matching our branching ratio to LEP data yields an 
octet matrix element with substantially smaller theoretical uncertainties than from hadronic 
processes |3|-|6| . The differential decay rate is dramatically softer than previous calculations. 
The small-;? peak in the differential distribution represents a clear signature of the octet 
mechanism that we regard as strong motivation for continued analyses by the LEP experi- 
mental groups. A measurement of the first moment would be a particularly interesting. 
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